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Abstract 
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tensor product graph methods which were originally formulated for quantum 
algebras. Connections with the Fateev-Zamolodchikov model are discussed. 
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1 Introduction 

The Yang-Baxter equation in its spectral parameter-dependent form has well-known 
applications in a variety of areas of mathematical physics, including exactly solvable 
classical models of two-dimensional statistical mechanics and integrable quantum sys- 
tems. The origins of the field trace back to the influential works by McGuire con- 
cerning many-body systems with delta function interactions [1], Yang's solution for 
interacting fermions [2], and Baxter's solution of the eight- vertex model [3]. Since 
then numerous solutions have been obtained including those leading to the Perk- 
Schultz model [4], the Andrews-Baxter- Forrester model [5], the chiral Potts model 
[6, 7] and the Hubbard model [8, 9] to name a few. 

One systematic method for solving the Yang-Baxter equation was found with the 
advent of affine quantum algebras, which are deformations of the universal envelop- 
ing algebras of affine classical Lie algebras [10, 11]. Through the Drinfeld double 
construction [10], the affine quantum algebras can be seen to belong to a larger class 
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of algebras which are referred to as quasi-triangular Hopf algebras. The double con- 
struction provides a universal prescription for constructing a quasi-triangular Hopf 
algebra, denoted D(H), from any Hopf algebra H and its dual algebra H*. Via this 
construction there exists a canonical element in the tensor product algebra of the 
Drinfeld double, R G D(H) ® D(H), known as the universal i?-matrix. This element 
provides an algebraic solution of the Yang-Baxter equation. Thus for each matrix rep- 
resentation of D(H), a matrix solution of the Yang-Baxter equation is automatically 
obtained. Solutions of the Yang-Baxter equation with spectral parameter arise when- 
ever matrix representations of D(H) can be parameterised by one or more continuous 
variables. In the context of affine quantum algebras, the spectral parameter arises 
naturally by consideration of an evaluation homomorphism from the affine quantum 
algebra to its non- affine counterpart. In this manner the representations of the affine 
quantum algebras which give solutions with spectral parameter are the loop repre- 
sentations. In principle the evaluation homomorphism provides a route to algebraic 
solutions of the Yang-Baxter equation with spectral parameter. However it is a tech- 
nically challenging task to undertake and has only been implemented explicitly for 
some low rank cases [12, 13]. 

Often, a more tractable approach is to start with a constant matrix solution of 
the Yang-Baxter equation and then ask the question of whether a spectral parameter- 
dependent generalisation exists, a procedure which is commonly referred to as Bax- 
terisation. For example, in cases where the solution of the Yang-Baxter equation 
without spectral parameter affords a representation of the Hecke algebra, there exists 
a procedure for constructing the associated spectral parameter-dependent solution 
[14]. This also applies in cases where the constant solution gives rise to a represen- 
tation of the Birman-Wenzl-Murakami algebra [15]. For solutions associated with 
representations of quantum algebras, a general prescription for Baxterising constant 
solutions of the Yang-Baxter equation (subject to satisfying certain conditions) was 
developed in [16, 17, 18] using the notion of the tensor product graph method. 

An alternative route to building solutions of the Yang-Baxter equation is in the 
framework of descendants. The most notable instances were found by Bazhanov 
and Stroganov [19], who showed that the chiral Potts model is a descendant of the 
six- vertex model, and Hasegawa and Yamada [20], who showed that the Kashiwara- 
Miwa model is a descendant of the zero-field eight-vertex model. In both cases the 
connection was found by constructing an appropriate algebraic structure through an 
L-operator. In [19] the algebraic structure is closely related to that of the quantum al- 
gebra U q (sl(2)) with q a root of unity, and in [20] the L-operator is expressed in terms 
of the Sklyanin algebra [21]. In this approach use of the L-operator allows for the con- 
struction of descendants through higher dimensional representations of the algebra. 
It is known [22] that the six-vertex and the zero-field eight-vertex models intersect at 
the zero-field six-vertex model whose descendant is the Fateev-Zamolodchikov model 
[23], which is precisely the intersection between the chiral Potts and Kashiwara-Miwa 
models. 

Another principal class of quasi-triangular Hopf algebras is the class of Drinfeld 
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doubles of finite group algebras. These algebras have recently received attention in 
relation to the description of anyonic symmetries in quantum systems, specifically in 
terms of the braiding properties of anyonic quasiparticle excitations which are de- 
scribed by constant solutions of the Yang-Baxter equation [24, 25, 26, 27]. However 
these algebras these have not been investigated in the context of solving the Yang- 
Baxter equation with spectral parameter, apart from some preliminary investigations 
reported in [28, 29]. Arguably the simplest family of Drinfeld double algebras for finite 
groups is that associated with the dihedral groups, denoted D n . The Drinfeld doubles 
of these algebras, denoted D(D n ), are finite-dimensional quasi-triangular Hopf alge- 
bras with a finite number of irreducible representations. For D(D n ) all the irreducible 
representations are known [28]. In the case where n is odd, the irreducible represen- 
tations can only have dimensions 1, 2 or n, whereas for n even the dimensions are 1, 
2 or n/2. The constant i?-matrices associated with the 2-dimensional representations 
can be Baxterised to reproduce the six- vertex model at roots of unity [28]. 

Our goal below is to obtain spectral parameter-dependent solutions of the Yang- 
Baxter equation associated with the n-dimensional (n odd) and n/2-dimensional (n 
even) representations. The general strategy we follow is to view these solutions as de- 
scendants of the six-vertex model, by firstly determining the appropriate L-operators. 
Once this is achieved, the technical aspects for determining the explicit form of the de- 
scendants is accommodated by adapting the tensor product graph method developed 
for quantum algebras [16, 17, 18] to the present case. It turns out that solutions we 
obtain for n odd are limiting cases of the Fateev-Zamolodchikov models [30]. However 
we stress that the anticipated U q (sl(2)) symmetry with q a root of unity degenerates 
in this limit, with the D(D n ) symmetry emerging. 

2 Preliminaries 

We first define some notation which is used extensively throughout this article. We 
often work modulo n G N, and hence define the following map: given a G Z then 
< a < n — lis the integer which satisfies a = a (mod n). We also use the following 
two delta functions: 



We use eij to denote an elementary matrix in M nxn (C) whose indices are considered 
modulo n. We adopt the convention that eo,o corresponds to the matrix with an entry 
in the nth row and nth column. These matrices obey the relation 




1, i = j, 
0, i±h 




1, i = j (mod n), 
0, i ^ j (mod n). 



We use the following convention for the product symbol: 
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2.1 The Yang-Baxter equation with spectral parameter 

We will use three variants of the Yang-Baxter equation (YBE), which is a non-linear 
matrix equation in End (V <8> V <E> V) for some vector space V. The first of these forms 
we call the constant YBE, given by 

R12R13R23 — -R23-R13-R12, 

where the subscripts of R refer to which vector spaces the operator is acting upon. 
That is, given 

R = ^ a% ® bi G End (V <g> V), 

i 

we have 

R12 = ^2 a i ® &i <8> -R13 = (g> / (g) 6j, etc. 

i i 

Here / denotes the identity matrix. The second form is known as the spectral 
parameter-dependent YBE and is given by 

R 12 {x)R 13 {xy)R 23 (y) = R 23 (y)R 13 (xy)R 12 (x), (1) 

where x,y G C. Given a parameter-dependent solution, R(z), constant solutions are 
recovered when we take the limits z — > 0, 1 and 00. We refer to invertible solutions 
of either form as R-matrices. As this article deals primarily with the parameter- 
dependent YBE we shall henceforth refer to it simply as the Yang-Baxter equation. 

Now suppose r(z) G End (V <g> V) and R(z) G End (W ® W) are i?-matrices with 
dim V < dimly. We describe R(z) as a descendant of r(z) provided that there exists 
a non-trivial invertible operator L(z) G End (V <8> W) (referred to as an L-operator) 
which satisfies 

r\2{xy' 1 )L l3 (x)L 23 {y) = L 23 {y)L l3 (x)r l 2{xy~ 1 ) ) Wx,y G C, (2) 

and 

Li2{x)L n (y)R 23 {x~ 1 y) = R23(x~ 1 y)L 13 (y)L 12 (x), Vx,y G C. (3) 

In this article we start with an i?-matrix associated with the zero-field six-vertex 
model and use the Drinfeld doubles of dihedral groups to generate L-operators and 
descendants. 

The third variant of the YBE used in this paper is 

Ri2(x)R 23 (xy)R 12 (y) = R2z{y)Ri2{xy)R,2z{x). 

Given a matrix solution R(z) to Equation (1), a solution to this form of the YBE is 
given by 

R(z) =PR(z), 
where P is the permutation operator, i.e. 

P = ^ e it j <8> e^j. 
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2.2 Drinfeld doubles of finite groups 



Next we recall relevant results from the Drinfeld double construction applied to finite 
group algebras, following [31, 32]. Given a group G, with identity e, we consider the 
algebra CG, whose basis vectors are the elements of the group. The multiplication 
and unit of CG are inherited from the group in the natural way. We equip CG with 
a coproduct, counit and antipode defined respectively by 

A(g)=g®g, e(0) = l and 7 (» = g~\ VgeG. 

With these maps CG becomes a Hopf algebra. We next consider the dual space of 
CG, 

(CG) = C{g*\geG}. 
The multiplication and unit are given, respectively, by 

g*h* = 5 h g and «(1) 

The costructure and antipode are defined by 
A(g*) = J2(hgT ® (O*, e(g*) = S e g and 

heG 

Under these maps (CG) Q is also a Hopf algebra. Using the dual and the original 
algebra we can construct the Drinfeld double, 

D(G) = C{gh*\g,he G}. 

We impose the relation 

h*g = gig-'hg)* 

and adopt the required remaining structure from CG and (CG) . Furthermore it is 
known that D{G) is a quasi-triangular Hopf algebra, containing the canonical element 

geG 

This element satisfies the following relations: 

KA(a) = A T (a)n, Ma e D(G), 

(id®A)ft = K 13 K 12 , 

where A T denotes the opposite coproduct. Consequently matrix representations of 
1Z provide solutions of the constant Yang-Baxter equation. 

In this paper we use the dihedral group D n , which is the symmetry group of a 
regular polygon with n vertices. That is, 

D n = {a, r\a n = r 2 = e, ara = t}, 

and has order In. For ease of calculation we divide the description of the representa- 
tion theory into different cases. 



g&G 

1(9*) = (g- 1 )*, Vg e G. 
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2.3 Representations of D(D n ) 
2.3.1 The case when n is odd 

We first consider D(D n ) for the case where n is odd. As stated in [28] the represen- 
tations of the double of a group are naturally partitioned by the conjugacy classes of 
the group. For these representations we consider w G C to be a primitive nth root of 
unity. Then the irreducible representations (irreps) are given in Table 1 below: 



Irrep tt 


Constraints 


7r((T) 


TT(r) 


ir(g*), geD n 






1 








1 < k < ^ 


fw k \ 
\ w~ k ) 


(":) 


( S » °) 
V° ^ 




1 <l < 
< k < n- 1 


fw k \ 
V w~ k ) 


CD 








2-a=l e i+l,i 




d g e i+i,i+i 



Table 1: The irreps of D(D n ) when n is odd. 

Here irf,^'^ and ir^ have dimensions 1, 2 and n respectively, and their associated 
modules are denoted Vf, V 2 (a ' b) and V±. 

Also associated with these representations are the solutions of the Yang-Baxter 
equation arising from the canonical element. If we apply a two-dimensional irrep we 
find the i?-matrix 



4> k) )n 



( 


w kl 
















w -kl 
















— hi 

w 





\ 














It was shown in [28] that this constant solution leads to the parameter-dependent 
solution 



/ w 2kl z 1 — ir 



-2kl. 



r(z) 



\ 











z^ 1 — z 






z^ 1 — z 










\ 



(4) 



w 2kl z 1 — w 2kl z ) 



This i?-matrix corresponds to the six-vertex model with zero-field. Also of interest is 
the representation of the canonical element using the n-dimensional irreps, 



n-1 



fan & ;r 



i,j=0 



In our construction of the descendant of r(z) we will use a linear combination of 
projection operators. It is known [32] that there exist operators which project D{D n ) 
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onto its ideals. These projection operators are defined by 

E « = w E Xa^g^gh*, 
a,heG 

where a is an irrep, d[a] its dimension and Xa is the group character defined by 

X a (a) = tr 7r Q (a), Va G £>(£>„). 

We also consider the projection operators 

p« = (^®^)A(S a ). (5) 

These operators project from 7r^ ® 7r^ onto copies of the irrep associated with a in 
the decomposition of tt^ (8) 7r^. Note that we consider 7r+ (g) 7r+ and ir~ (g) 7r~ together 
since 

(tt+ ® vr+)A(a) = (n- ® Tr-)A(a), Va G D(D n ). 
We calculate p° explicitly, using the expression given in Equation (5), obtaining 

Pn = E E X«((^') V 1 )^^)^-,- ® ^(<?Ki- 

g<=D n i,j=0 

From this we can see that the only irreps with non-zero projection operators will be 
associated with the conjugacy classes {e} and {a 1 , cr~ 1 } for 1 < i < r± ^- This implies 
that Cg> 7r^ decomposes only into one and two-dimensional irreps. 

For convenience, we slightly modify our notation for the irreps. Instead of using 
and Ti2 ,k \ we use ordered pairs corresponding only to irreps that appear in the 
direct sum decomposition of (g) ir^- The correspondence is summarised in Table 2: 



a 


Irrep 


Constraint on a 


Constraint on 6 


(0,0) 








(0,6) 






1 < 6 < L^J 


(0,6) 


(0,n-2&) 
^2 




[n±3j < & < n=l 


(a, 6) 


(2a,26) 
^2 


1 < « < L¥J 


< 6 < n - 1 


(a, 6) 


(n-2a,n-2b) 
^2 




< 6 < n - 1 



Table 2: The ordered pairs labelling the irreps for D(D n ), n odd. 

Here [a\ denotes the floor of a. Then the projection operator for irrep a = (a, 6) is 
given by 

Pn = — E [w^i+a+W+o ® e *+W + W 2k 'ei-a+jj-a ® Cj+^i] 
77/ 

where 
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c a = \h « = (0,0), 
\l, a ^(0,0). 

We have calculated non-zero projection operators for (n 2 — l)/2 two-dimensional 
irreps and one 1-dimensional irrep. Using a counting argument, it is clear that these 
provide a complete decomposition of <g> 7r^ into irreps for n odd. 

2.3.2 The case when n is even 

We similarly catalogue the irreps of D(D n ) for the case when n is even. We set 
n = 2m and let w G C be a primitive 2mth root of unity. 



Irrep 7r 


Constraints 


7r(CT) 


7T(r) 


7T(#*), ^ G £> 2 m 


"l,e 


a,&G {0,1} 




(-1)- 




71 l,<7 m 


a,be {0,1} 










1 < k < m 


/w fc \ 
\ w- k ) 


C o) 


V° ^ 




1 < k < m 


\ w~ k ) 




Vo <7 


(l,k) 


1 < / < m - 1, 
< fc < 2m - 1 


fw k \ 
\ w~ k ) 




Vo ^7 




a,be {0,1} 


Er=i(-i) a5 ^v-i 






" m,(TT 


a,&G {0,1} 


Er=i(-i) a5 ^M-i 


(-i) 6 Er =1 ^ 


g £k+l,k+l 



Table 3: The irreps of D(D 2m ). 



We again concern ourselves with the two-dimensional representations applied to 
the canonical element. We find the constant solution 



( 


w kl 


















w -kl 


















— hi 

w 







V 











w kl , 


) 



which leads to the parameter-dependent solution 

/ w 2kl z~ x - w- 2kl z 

( \ 
r{z)= Q 

\ 



z^ 1 — z 

o 



o 






w 2kl z " — w~ 



l z ) 
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This result only differs to that of D(D n ) where n is odd by the possible choices of the 
root of unity. Note also the similarity in the representation of the canonical element, 
given by 

TO — 1 

(7T m <g> ir m )n = ^ e i+j,i~j ® e i,ii 

i,j=0 

when 7r m is either TVm'P or Tim^X for b G {0, 1}. 

We again use Equation (5) to derive our projection operators. We first look at 
the case when m = n/2 is odd. We find that 

jr i _ m-1 
g&G i,j=0 

where 7i m is any one of the m-dimensional irreps. As in the case of D(D n ) for odd 
n, we introduce an ordered pair notation for the irreps that appear in the direct sum 
decomposition of KmP <8> 7rm'r and TTm$ T ® 77™,'^ for be {0, 1}. The correspondence 
is summarised in Table 4: 



a 


Irrep 


Constraint on a 


Constraint on 6 


(0,0) 
(0,6) 
(a, 6) 


-(0,26) 
^2 

(2a,26) 
"2 


1< a < m " 1 


1< 6 < m " 1 
< b < m - 1 



Table 4: The ordered pairs labelling the irreps for D(D 2m ), m odd. 



Then the projection operator associated with irrep a = (a, b) is given by 
where 



m 



c a = ^ 2 ' 



|, a = (0,0), 
1, a ^(0,0). 



We note that these projection operators here match those arising in the case of 
D(D n ) where is n odd. 

Similarly, when m = n/2 is even we have the projection operators 

Pm = Xa(((r 2j )*g~ 1 )nm(g)e i - j , i - j <g> 7T m (#)e M , 

<?eG i,j=0 

for any m-dimensional irrep 7r m . The ordered pair notation for the irreps occuring in 
the decomposition of TTm'fr <8> TTm'Jr for b G {0, 1} is given in Table 5, whereas those 
occuring in the decomposition of iTm'$ T ® vrm^l for 5 G {0, 1} are given in Table 6. 
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a 


Irrep 


Constraint on a 


Constraint on 6 


(0,6f) 
(f,6f) 

(0,6) 
(f,6) 
(a, 6) 


(0,6) 
"l.e 
(0,6) 

" l,cr m 

(0,26) 
^2 

(m,26) 
n 2 

(2a,26) 
"2 


1 < a < f - 1 


b €{0,1} 

&e {0,1} 
1 < b < f - 1 

1 < b < f - 1 
< 6 < m - 1 



Table 5: The ordered pairs labelling the irreps for Z)(D 2 m), m even, case (i). 



a 


Irrep 


Constraint on a 


Constraint on 6 


(0,6?) 
(f,6f) 
(0,6) 

(f,6) 
(a, 6) 


"l.e 

" l,cr m 

(0,26) 

(m,2b) 

(2o,26) 
"2 


1 < a < f - 1 


6 €{0,1} 
6 G{0,1} 
1<6< f -1 
1 < 6 < f - 1 
< 6 < m - 1 



Table 6: The ordered pairs labelling the irreps for D(D 2m ), m even, case (ii). 



We obtain the same projection operators irrespective of which of the two cases we 
consider. In particular, the projection operator corresponding to the irrep a = (a, 6) 
is 



m-l 

C 



m 

i, 3 =Q 



where 



1, a G {0, f } and 6 G {0, f }, 



c 



otherwise. 



Again, a counting argument verifies that these are the only non-zero projection 
operators. As the projection operators are the same for the different m-dimensional 
irreps, we can without loss of generality consider only -Km'P- 



3 Descendants associated with D(D n ) - n odd 

In this section we construct a family of solutions of the Yang-Baxter equation using 
D(D n ), where n is odd and n > 2. 



3.1 Construction of the L-operator 

We begin by constructing an operator L(z) G End(V 2 ^' fc ' ) <8> V^)- We observe that 
V^ 1 '^ <8> = © V~, suggesting that the L-operator will be a linear combination 
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of two intertwining operators with one parameter-dependent coefficient (after scaling). 
Hence we adopt the following ansatz for the L-operator: 

L(z) = (4 k) ®ir±)[n + h(z)(n T y 1 ] 

where (JZ T )~ l = ^g* <S> g -1 - Applying the representations leads to the following 

g&G 

expression for L(z): 

n-1 

L(z) = {{w 2{i ~ 1)k ei,i + w- 2( - i ~ 1)k e 2:1 ) <g> e iti + h(z) [e ltl ® e*_ M + e 2j2 ® e i+M ]} . 

i=0 

Performing a basis transformation which leaves r(z) invariant, we obtain 

n-1 

L(z) = {{w 2ik e l)2 + w~ 2ik e 2 ,i) ® e M + h(z) [e 1A ® + e 2 , 2 ® e i+J)i ]} . 

i=0 

Substituting r(z) and L(z) into Equation (2), we find only one constraint on h(z), 
namely 

h(x)y - /i(j/)a; = 0, 

which has the solution 

h{z) = Cz, 

for any C G C. We are free to rescale the parameter z without affecting the descen- 
dants given in the next section, so without loss of generality we can choose C — 1. 
We have therefore shown the following: 

Proposition 3.1. The L-operator given explicitly by 

n—l 

L(z) = Y { (w 2ik e h2 + w- 2ik e 2 ,i) ® e M + z [e 1A ® + e 2>2 ® e i+lti \ } , (6) 

i=0 

and the r(z) given in Equation (4) together satisfy Equation (2). 

It is important to comment that, up to a basis transformation, the L-operator (6) 
is a particular limit of the general L-operator discussed in [19] in relation to the chiral 
Potts model. However, the associated U q (sl(2)) structure with q 3 = 1 described in 
[19] is lost in this limit. Specifically, Equation (2.9) in [19] does not hold. Amongst 
the defining relations of the generalised U q (sl(2)) algebra this means, in particular, 
that 

[e, /] = 

where e and / are the raising and lowering generators respectively. The above relation 
is indicative of the fact that the U q (sl(2)) structure degenerates in this limit. Instead 
we find that the symmetry algebra of D(D n ) emerges when the L-operator is given 
by Equation (6). 
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3.2 Construction of the descendants 

Here we put forth a predicted form of a descendant, R(z), and determine the con- 
straints on this form. We then impose additional constraints enforcing that the de- 
scendant inherits certain properties from representations of the canonical element 
from D(D n ). For convenience we shall use an alternate form of Equation (3): 

J R 2 3(^- 1 )L 13 (x)L 12 (y) = L 13 (y)L 12 (x)R 23 (xy~ 1 ), (7) 

where 

R(z) = PR(z) 

and L(z) is given by Equation (6). As the descendants R(z) commute with the action 
of the coproduct, we assume they are of the form 

R(z) = J2Mz)p a , (8) 

where p a are the projection operators previously calculated, f a (z) are continuous 
functions and S is the set of ordered pairs which correspond to non-zero projection 
operators. For convenience we rescale our projection operators, henceforth using 

n-1 

P a = ^2 [ w23be i+a+i,i+a ® e i+j>i + w~ 2h ' J ' ei_ a+j ^ a ® e i+j4 ], a = (a, b) e S. 

i,j=0 

Applying p 23 and p 23 to the left- and right-hand sides respectively of Equation (7), 
we find 

f a (xy~ 1 )p 23 L 13 (x)L 12 (y)p^ 3 = f l3 (xy' 1 )p^ 3 L l3 (y)L l2 (x)p^ 3 . (9) 
We calculate that 

L 13 (x)L 12 (y) = en <g> A(x, y) + e 12 ® B(x, y) + e 21 <g> C(x, y) + e 22 <8> D(x, y), 
where 

n-1 

M x > V) = XI w2b ~ t)ke i,i ® ej-j + xy e;_ M <g> e^ hj , 

i,j=0 
n-1 

B(x, y) = ^ :ru> 2lfc e M <g> e^-jj + yw 2jk e i+ i ti ® e^-, 

n-1 
n-1 

^(z, 2/) = X w2(l_j)fee i,i ® e i,i + ^ e *+M ® e i+/,r ( 10 ) 

i,j=0 
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Substituting these functions into Equation (9), we find the following constraints on 



f a {xy 1 )p a A(x,y)p^ 
f a (xy~ l )p a B{x,y)f 
f Q (xy- 1 )p a C(x,y)p^ 
fa{xy~ 1 )p a D(x,y)p^ 



fp{xy 1 )p a A(x,y)p p , 
fp(xy~ 1 )p a B( y y,x)p 13 , 
fp{xy- l )p a C{y,x)p^ 
f^(xy~ 1 )p a D(x,y)p IB . 



(11) 
(12) 

(13) 

(14) 



Note that Equations (11,14) are automatically satisfied. 

The above relations provide a tensor product graph [16, 17, 18]. Assigning each 
irrep to a vertex of a graph, we connect the vertices labelled a and f3 by an edge if 

P a x(x,y)p p ^ 

for either \ = B or x = C. The tensor product graph for the case I = k = 1 is 
depicted in Figure 1. An edge connecting vertices a and /3 signifies that the functions 
f a (x) and f/3(x) are constrained by Equations (12,13). 



(0,0) 



(0,1) 





(l,a±i) 



i i — 4 



2 ' 2 ' ^ 2 ' 2 ^ 2 ' 2 



i i 



(2=A,0) (3=1, n _i) 



Figure 1: Tensor product graph for odd n when I — k — 1. 

We now give a series of propositions and lemmas which result in a solution to 
Equation (7). 

Proposition 3.2. For R(z) defined by Equation (8) to be a solution of Equation (7), 
the following four constraints must be satisfied: 

= 5 d k J C a + l [ha, b) {z){ Z W^- 1 -^ + 1) - f (c4) (z)(™ 2i( - l - a)k - bl) + *)] , 

o = s d k - b s l - a [f ia M^ 2{i ' l+a)k+bl) + 1) - f(cM^ 2{{ - l+a)k+bl) + Z)] , 

= 5^ b 5 l + c a {f {a>h) {z){zw 2 ^-^-W + 1) - f (c4) (z)(w 2 «- l - a ^ + Z)] , 

= [f(a, b) (^ 2{i ' l+a)k+bl) + 1) - hc4m^ 2{( - l+a)k+U) + Z)] , 

where the pairs (a, 6) and (c, d) belong to S. 
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Proof. We consider constraint Equation (12). Given the operators 

n-l 

B (x,y) = ^ xw2lke i,i ® e j-i,j + yw 2]k e i+lti <g> e jd , 

i,j=0 

and 

n-l 

P a =^2 l w2jbe i+a+j,i+a ® e i+jti + w~ 2b3 ei- a+jti - a <g> e i+ ^], 

i,j=0 

we first calculate 

p {a ' b) B{x,y) 

n-l 

= ^ [w~ 2b3 e^ a+jti _ a <g> ei + j,< + w 2fej e i+a+i> j +a <g> e i+j>i ] [xw 2sk e StS <g> e t _ M + yw 2th e s+ i jS <g> e t;t ] 

j,jr',S,t = 

n-l 

= £ [x«; 2 «*-'- a ) fc - 6 ^e t _^ a+J - t _ z _ a ® e t _ i+J - 1 + x W 2 «'-' +a ) fc+ ^e^ +a+J , t _, +a ® e t _ J+j) t 

+ yw 2 (* fe - 6 ^e t _ a+jit _ a _i <g> et+j-t + 2/u> 2( ' fc+6j) e t +a+^+a-« <g> e t+J - t ]. 
This leads to 

p^B(x,y)p^ d) 

n-l 

= J2 W at ^ a)k ~ b3) et-i-a + ,,t-i-a ® e t - l+jtt + xw 2 ^- l+a ^e t . l+a+ ^ l+a ® e t _, +j , t 

+yw m-bj) et _ a+j t _ a _ l et+ . t + ^ 2 ( tfc+ ^ e , +a+J - t+a _, ® e t+J - 

n-l 

= [«i; 2 « t - , - a) ^'-* ) ^e t _ I _^- l4 _ c _ 1 , ® e t _^, t _„ 

j,t,u=0 

+ a:u ,2((*-'+-)^-*)5'-« et _ J+o+ et _ J+J . t _„ 

+ j /u ,2( tfc -6j-d V ) 5 a+i et _ a+ ^ .^^ + yw m+bj-dv)^-a et+a+jt _ c _ v ^ ^ . ^ 

+xw 2 ^ t - l ' a ^ +d ^5 l X a e t - l - a+3 , t+c - v ® e t _ J+ji t _„ 

+ XW 2(( t -i + o)fc + 6j+d V ) 5 a-i et _ j+o+ et _ J+ . )t _ w 

+J/w 2( tfc -6 J - + d V ) 5 a+i et _ a+ .^^ + y w m+bj + dv)^a-l et+a+jt+c _ v ^ ^ . 

n—l s 

= nJ2{ ™ 2{dV ~ M) [xw 2 ^ 1 -^ + y] ¥ k+b 5 c a+l e 3 ^ c ® e,, 
+ w 2( ^ } ^((-'^WO + y ] ¥ k _ b 8 l re j+a , v - c ® 

+w -(2*,+6i) [^((-J-aJk-W) + y ] J-^^e ._ a ^ +c (g, e . v 
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and finally arrive at the desired four constraints. 

Very similar calculations for Equation (13) leads to the same four constraint equa- 
tions. 



Lemma 3.3. Let w 2 be a primitive nth root of unity and I, k be integers such that 
gcd(l,n) = gcd(k,n) = 1. Then 



for all a G N, b G Z. 

The proof is omitted. Henceforth we choose I, k satisfying gcd(l, n) = gcd(k, n) = 1 
so that we may use the above lemma. 

Proposition 3.4. Let S' = {(a, b)\a > 0,6 G Z} ; and consider a set of functions 
{f(a,b)\{ a ,b) G S'}. If the functions satisfy the relations 



then the functions also satisfy the constraints given in Proposition 3.2 for all (a, b) G 
S. Moreover, every set of functions satisfying the conditions of Proposition 3.2 can 
be extended in a unique way to a set of functions defined on S' satisfying the above 
conditions. Hence the two sets of constraints are equivalent. 

The proof is straightforward and omitted. 

Lemma 3.5. The set of functions 



□ 





f(a,b)( Z ) 

f(a,b)( Z ) 
f{a,b)(z) 




f( n -a,-b)(z), (a,b)eS', 




satisfies the conditions in Proposition 3.4 given 



f(0,b)( Z ) = f(0,-b)( z ) = /(0,6+n), 



for b G Z. 
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The proof follows from Lemma 3.3. 

Substituting the functions given in Lemma 3.5 and their associated projection 
operators into our form for R(z) gives the operator 



n-l 



al 



*M = n E ^ n I 



i,j,a,b=0 



p=l 



_|_ UJ 2/((2p-l)fe+6-aM- 1 ) 
+ ^2/((2p-l)fc+6-afc/- 1 ) 



/(0,fe-afc«- 1 )(^) e i+a+j,i+a ® e i+i,i 5 



which satisfies Equation (7). Note that in these calculations the only property that 
we have used is that w 2 is a primitive nth root of unity, irrespective of whether n is 
odd or even. This allows us to use these calculations later on without alteration. 

For R(z) to be a descendent, it must also satisfy the Yang-Baxter equation. For 
the moment we will set I — k — 1, and consider the more general case in Section 5.1. 
Then the functions become 



f ( \_fj ( z + wm ~ 1+b) \ / < \ 

}(a,a+b){Z) - [I y - ~ znj2(2] _ l+b) J 7(0,6) (Z). 



3=1 



(15) 



This gives the operator 



n-l 
2 a 



R(z) 



a=0 j=l 



z + w 2{2 ^ 



p 



(a,o) 



n-l 
2 



6=1 



a=l 



n 



z + w 2( 2j -l+6) 
1 + ZW 2(2j-l+6) 



P 



(a,a+6) 



„2(2j-l-6) 



+ zw m-i-b) 



p 



(a, a— 6) 



or equivalently 



n— 1 



i,j,a=0 



n— 1 a /• 

6=0 p=l v 



z _|_ ^,2(2^-1+6-0) 



+ 2U) 2(2p-l+6-a) 



f(o,b-a)(z) 



&i+a+j,i+a® e i+j,i- (16) 



In there remain arbitrary functions, one for each disconnected component 
in the tensor product graph of Figure 1. This is in some contrast to the examples 
studied in [16, 17, 18] for which the tensor product graphs are connected. To ensure 
R(z) inherits properties from the representation of the canonical element, we enforce 
additional conditions. Firstly, we impose that R(z) is self-adjoint, which is equivalent 
to 

R 12 (z) = Rl^z), zeR, 
where f is the adjoint operator. 

Proposition 3.6. The following three statements about the matrix R(z) of Equation 
(16) and coefficient functions of Equation (15) are equivalent: 
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(i) R l2 (z)=R\ l (z). 

(ii) (f {a , b) (z))* = f {a ,_ b) (z), V(a,b)eS>. 

(Hi) (f{o,b)(z))* = f(o,b)( z ) = f(o,b+2c)(z), Vz E R and b,ceZ. 
Here * denotes complex conjugation. 

Proof. Assume (i) holds. Looking at the general form of the projection operators 
we find that 

blfV = (pSt 6) )- 

Recalling the form of R(z) given in Equation (8) and using the linear independence 
of the projection operators, this implies that 

(/(«*)(*))* = /<•,-*)(*), V(a,b)eS'. 

Conversely, (i) follows directly from (ii), and hence statements (i) and (ii) are 
equivalent. 

Now assume (ii) holds. Equation (15) gives us 



and 



We see that 



t(i,i+b){z) - y l + zw 2(b+i) ) no,b){ z ), 

}(l,-l-b)(z) = y i + zw -2(b+l) ) J(0,-b-2)(z). 



/ z _|_ w 2(6+l) \ / z + ^2(6+1) \ 

(/(1,-1-b)^))* = (^ 1 + ^2(iH-l) J (f(0,-b-2)(z))* = y i + zw 2(b+l) J f(0,b+2)(z), 



which implies that 

f(o,b)(z) = f(0,b+2)(z). 
Combining this with the constraint of /( ,&) (z) we find 

(/(0,6)(*))* = f(o,b)( z ) = f(o,b+2c)(z), 

for b,c E Z. 

Conversely, suppose (Hi) holds. It follows that 



f m - TT + ^ 2(2j " 1+b ' a) A f M 

J(a,6)W - 11 ^ 1 + ^2(2j-l+ft-a) J /(0,6-o)^J- 
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Hence 



/ " / z + w 2(2j-l+b-a) \ 
(/(<»,&)(*))* = ^ + zw 2(2j-l+b-a) ) f(P,b-a)(z) 

= II ( 1 + zw -2(2j-l+b-a) J (/(0,&-«)(*))* 

" / z _|_ w 2(2(-j)+l-6+a) \ 
= 11 ^ 1+ ^2(2(-i)+l-N-a) J /(0,0-6-a)^) 

11 l^ 1+ ^2(2,-l-6-a) J t(0,-b-a){z) 



3=1 
f(a-b)(z). 



This completes the proof. 



□ 



As the second index of the function can be considered modulo n and n is odd, im- 
posing self-adjointness implies that there is only one arbitrary function left. Without 
loss of generality we consider it to be /(o,o)(z). This can be seen as an overall scalar 
of our operator thus we are able to set it to a constant. Now we impose the limiting 
condition 

WmR{z) = ±{^®^)TZ. (17) 
This sets the scalar to be f(ofl)(z) = 1, so the functions become 

f(a,b){ Z ) - I M l + zw2{2j _ 1+b _ a) I • 
j=l \ / 

Hence 



*] = Ef"'' )+ EEn i 

6=0 a=l 6=0 i=l 

or equivalently 



Vn-l a /z + w2(2j _ 1+b _ a) 



+ ZW 2(2i-l+6-a) 



n-1 



i,j,a=0 



1 n_1 a / 7 _l ,„2(2p-l+6-o) \ 

11 I 1 + zw 2(2p-l+b-a) I 
h—n r>=1 V / 



&i+a+j,i+a ® 



6=0 p=l 

To simplify this we define the functions 

1 «-! i a / , 2(2p-l+6-o) \ 

«*>(*> - i £^/ ( »,*,w = i n ■ 

6=0 6=0 p=l v 
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where a G N and j G Z. In terms of these new functions 9( a ,j)( z )j our operator 
becomes 

n-1 

R( z ) = ^ 9{a,j)( z )^i+a+j,i+a ® 
i,j,a=0 

or 

n-1 

R ( z ) = Yl 9(a,j)(z)e i+j ,i +a ® e i+a+j)i . (18) 

i,j,a=0 

Remark. The operator R(z) of Equation (18) satisfies the following properties: 

(i) R{z)* = R(z), Vz G R, 

(ii) R\z) = R(z), Mz G C, 

(iii) R-^z) = R(z), \/z G C, 

(iv) R 12 {z)R 21 {z- 1 ) = I <g> I, Vz G C, 

(v) limi?(z) = ±(tt± ® 7r±)ft, 

(vi) limi?(z) = P. 

z— >1 

Proposition 3.7. The operator R(z) of Equation (18) is a descendant if and only if 

n—1 n—1 

9(a,k-d) (x)g (k , a - b) (xy)g (b ^ k) (y) = ^ 9{c,k-b) (x)9(k,c-d) {xy)g^ a - k) (y) , 

k=0 k=0 

for < a, b, c, d < n — 1 . 

This follows directly from the Yang-Baxter equation. 

The matrix R(z) of Equation (18) agrees with the i?-matrix obtained in [28] when 
n — 3. It has also been verified computationally for all odd n < 17 that the functions 
g(z) satisfy the equation given in Proposition 3.7 above, and thus that R(z) is a 
solution to the Yang-Baxter equation. The computations were performed with z 
treated as an arbitrary complex number and w as an arbitrary primitive nth root of 
unity. 

Conjecture 3.8. The matrix R(z) given in Equation (18) is a descendant of the 
zero-field six- vertex model with D(D n ) symmetry for all odd n. 
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4 Descendants associated with D(D2 m ) 



In this section we construct a family of solutions of the Yang-Baxter equation using 
D(D 2m )i treating the cases where m is even and odd separately. It will be shown 
that although in both cases the representations have the same general form, the case 
where m is odd has more in common with the solutions from the previous section 
associated with D(D n ) where n is odd. 

Throughout this entire section, we consider w to be a primitive 2mth root of unity 
and we use the r(z) stated previously and associated with n 2 = tt^'^ ■ The L-operator 
arising from the cases of m odd or even can be considered at once, since the approach 
is similar to that of Section 3.1. 

4.1 Construction of the L-operator 

In this subsection we construct an operator L(z) G End(V 2 <S> V m ) using a similar 
approach as before in Section 3.1. Again, we assume that the L-operator is of the 
form 

L(z) = (n 2 (g ) 7r m )[R + h(z)(R T )- 1 ], 

where 7r m is either 7Tm,'? or iXm b X with b e {0, 1}. Applying this representation and 
an appropriate basis transformation on the two-dimensional space we find 

m—1 

L(z) = {(w 2ik ei,2 + w- 2ik e 2A ) <g> e M + h(z) [e M <g> e^ M + e 2 , 2 ® e i+M ]} . 

We note that the basis transformation is of the same general form used in Section 
3.1, as are the operators r(z) and L(z). Thus it follows that our L-operator is of the 
form 

m—1 

L(z) = J] {{w 2ik e 1)2 + w~ 2ik e 2A ) ® e i}i + z [e M ® e^ hi + e 2j2 ® e i+M ]} . 

4.2 Construction of the descendants when m is even 

Here we consider the case of D(D 2m ) where m is even. We use r(z) associated with 
with w a primitive 2mth root of unity. As before, we look for descendants of 
the form 

PR(z)=R(z) = Y l fMp a , 

where p a are the projection operators previously calculated, f a (z) are continuous 
functions and S is the set of ordered pairs which correspond to non-zero projection 
operators. We require that R(z) satisfies Equation (7). As before, we use rescaled 
projection operators: 

m 

P a = E [w 2:ib e i+a+jti+a (g) e i+jji + w- 2b3 ei_ a+j4 _ a ® e i+j4 ], a = (a, b) e S. 
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As these operators and the L-operator are equivalent to those found earlier, we can 
use the previous calculations to arrive at the tensor product graph shown in Figure 
2. 




Figure 2: Tensor product graph for even n when I = k = 1. 
We again define the set 

S' = {(a,b)\a e N,b G Z}, 

and extend the functions f a (z) to that set in a way analogous to that of Proposition 
3.4. We restrict ourselves to k,l satisfying gcd(7,m) =gcd(A;,m) = 1. Following the 
calculations of the previous section, we determine that if R(z) is a descendant then 
the functions f( a ,b)(z) must satisfy 

"Ll / z _|_ w 2l((2j-l)k+b) \ 

/(a.oW-i+ft)^) = 11 ( J T 2J((2j-l)fc+6) J fo,b)( Z )> 
3=1 V 7 

f(0,b)(z) = f(0,b+m)(z), 

f(p,b)(z) = f(o,-b)(z), V(a,b)eS'. 
We now choose to set I — k — 1; this gives the functions 

« / g + w 2 (2j --l + 6) \ 
/(a,«+6)(«; =H I j - ^2(2^-1+6) J /(0,6)W- 
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This produces the operator 



2 a / z + w 2(2 3 -l) 



a=0 j=l 

m 

2 a 



P 



(a, a) 



+/ ( o,f)W End 



+ u ,2(2,-l+f) 



21/' 



2(2j-i+f; 



a=0 j=l 

fm_ 1 



2 1 



6=1 



a=l 



n i 



2 + w 



2(2j-l-6) 



+ ^2(2^-1-6) 



P 



(a, a- 



+ n i 



z + w 2(2 J -l+f) 



P 



(a,a+f ) 



or equivalently 



m— 1 

E 

j,j',a=0 



m— 1 



6=0 



p=l 



2 + U 



,2(2p-l+6-a) 



+ ZU ,2(2p-l+6-a) 



/(0,6-a)(z) 



Here we again enforce that our operator is self-adjoint. We are able to use the previous 
calculations and recall that 

(7(0,6) W)* = f(0,b)( Z ) = f(0,b+2c)(z), 

\/z G K. and 6, c e N. The second index of the function can be considered modulo m 
and as m is even the functions are partitioned into two sets. We see that we have 
only two functions in which we have any freedom left; without loss of generality we 
consider them to be /(o,o)(<z) an d f(o,i)(z). As there are two functions we find that we 
cannot consider them an overall scalar, so we need to impose additional conditions. 
As the m-dimensional representations of the canonical element are self-adjoint, we 
enforce that R(z) is unitary . From this it follows that 

f(o,o)(z) = ±1 and / (0 ,i)(^) = ±1. 
Imposing the limiting condition given in Equation (17) sets 

/(o,o)(*) = /(o,i)(*) = 1- 

This yields the operator 

z + w 2(2p-l+6-a) 



m—1 



R(z) = E 

i,j,a=0 



1 m—1 a / 

6=0 p=l v 



zw 2(2p-l+b-a) 



^"i+a+j,i+a §9 



(19) 
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To simplify this we define the functions 



m—1 m—1 a 



1 _|_ zw 2(a-b+l-2k) 



) 



g { a,j)(z) = ^J2 w2bj M z ) = ^J2 w2b3 II 

b=0 b=0 k=l 



_|_ w 2(a-6+l-2fc) 



where a G N and j G Z. With the use of these new functions our operator becomes 



m—1 



(20) 



ij',o=0 



Remark. The operator -R(z) of Equation (20) satisfies the following properties: 



Note that the last property shows that R(z) does not satisfy regularity, unlike the 
R(z) constructed from D(D n ) when n is odd. We do, however, have the following 
analogue of Proposition 3.7: 

Proposition 4.1. The operator of Equation (20) is a descendant if and only if the 
following constraint is satisfied: 



for < a, b, c, d < m — 1 . 

It has been computationally verified that the functions g(z) satisfy the conditions 
in Proposition 4.1 above for all even m < 16, and hence that R(z) is a solution to 
the Yang-Baxter equation. 

Conjecture 4.2. The matrix R(z) given by Equation (20) is a descendent of the 
zero-field six- vertex model with D(D2 m ) symmetry for all even m. 



(i) R(z)* = R(z), Wz G R, 

(ii) R\z) = R(z), Vz G C, 

(iii) R-^z) = R(z), Vz G C, 



(iv) R 12 (z)R 21 (z- 1 ) = I <g> I, Vz G C, 

(iv) \\mR(z) = ±(tt± ® t±)H, 

(v) \mvR(z) ^ P. 



m—1 m—1 



X] 9(a,k-d) (x)g {k!a - b ) {xy)g (bfi - k ) (y) = ^ 9(c,k-b) (x)g {kjC - d) (xy)g (dA _ k) (y) , 
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4.3 Construction of the descendants when m is odd 

For completeness, we include the descendants associated with D(D 2 m) for odd m. In 
our construction of the descendant of r(z) we use a linear combination of projection 
operators found in Subsection 2.3.2. Using these projections, for a self-adjoint de- 
scendant of r(z) which limits to a representation of the canonical element, we obtain 
the operator 

m—l 

R(z) = £ 

i,j,a=0 

This is equivalent to the operator R(z) found using D(D m ) where m is odd. Hence 
it obeys the properties stated in the remark in Section 3.2. Moreover, R(z) is a 
descendant of the zero-field six-vertex model with D(D2 m ) symmetry for odd m < 17 
and conjectured to be a descendant for all odd m > 3. 



j m— 1 a / 

6=0 p=l v 



z _|_ w 2(2p-l+6-a) 



+ zw 2(2p-l+6-a) 



H+a+j,i- 



5 Generalised descendants 

In this section we investigate some of the choices we made in order to obtain de- 
scendants. Specifically, we investigate the effects of choosing a different initial 2- 
dimensional irrep and of imposing fewer conditions on the descendant. 



5.1 Dependency on choice of irreducible representations 

In earlier sections we constructed descendants by starting from the i?-matrix asso- 
ciated with the two-dimensional irrep tt^'^ ■ Here we explore the consequences of 
starting with a different two-dimensional irrep, namely any ir^ 1 '^ where gcd(Z,n) = 
gcd(/c,n) = 1. We consider n-dimensional descendants for any n > 3, with w a prim- 
itive nth root of unity. We recall that 

/ w kl z- 1 - w~ kl z \ 

, z^ 1 - z w kl - w~ kl 

r { z > ~ o w ki _ w -ki z -i _ z o 

\ w kl z- 1 - w~ kl z ) 

and 

n-l 

L ( z ) = £ {{w lk e^ 2 + w~ lk e 2 ,i) <8> e it i + z [e 1A ® e^ hi + e 2 , 2 ® e i+ i^} 

i=0 

where < I < | and < k < n — 1, although here w differs slightly from that used 
in Section 2.3. We now repeat our earlier calculations for our more general / and k, 
obtaining 

I n ~ l ( z _|_ w lk(2p-l+b-al- 1 ) \ 

k ^ = n S wh3k \{ [ l + zw ik { 2 P -i +b -ai-^) f(o,bk-aki-i)(z)e i+ a +j ,i +a ®e i+j , i . 

i,j,a,b=0 p=l \ J 
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We consider the change of basis on the n-dimensional space which yields 



-> ejj-i^-i, % E Z. 
Under this change of basis we find that our L-operator becomes 

n— 1 

L(z) = {(w ilk ei,2 + w~ ilk e 2 ,i) ® e M + z [e hl <g> + e 2>2 <g> e i+M ]} , 

i=0 

while -R(-z) becomes 

= - E ^ II [ 1 + zw ik(2p-i + b-a) ) f(o,( b -a)k)(z)e i+a+j:i+a ® e i+jii . 

i,j,a,b=0 p=l 

From this we see that these different choices of initial two-dimensional represen- 
tation provide equivalent descendants. The different choices yield different basis 
transformations, a permutation on the arbitrary functions and a change of the root 
of unity, which must remain primitive. Thus any 2-dimensional irrep satisfying 
gcd(/,n) = gcd(/c,n) = 1 results in equivalent descendants. This is unsurprising 
as R(z) is real whenever zeR, and choosing different /, k effectively just changes the 
root of unity being used. 



5.2 Descendants with an extra parameter associated with 
D(D2m) when m is even. 

Here we return to D(D 2m ) where m is even and we find more descendants by imposing 
fewer constraints. We use the general form 



m—1 



R(z) = E 

i ,j,a=0 



m—1 



rn 



6=0 



p=l 



_|_ w (2p-l+b-a) 

+ zw( 2 *'~ 1+b ~ a ') 



f(0,b-a)( 



&i+a+j,i+a ® 



where w is now a primitive mth root of unity. We recall that previously we imposed 
that R(z) is self-adjoint, unitary and obeys a limiting condition. If we ignore the 
limiting condition but enforce the other two conditions we find that 

/(o,26)(*) = f(o,o)(z) = ±1 and /( ,2&+i)(*) = f(o,i)(z) = ±1, 

for all b E Z. Without loss of generality we can set 

f m (z) = 1. 

This gives us two possibilities; the first choice is /(o,i)(z) = 1, which yields 



m—1 

R + {z) = E 

i,j,a=0 



m—1 



m 



b=0 



P =i 



_|_ ^Pp-l+ft-a) 
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(21) 



This operator corresponds to the conjectured descendant given in Equation 19. The 
second option is /(o,i)(z) = — 1, which gives the operator 



m— 1 
i,j,a=0 



m— 1 



rn 



£(-ir + v<n i 



f>=0 



P =i 



_|_ w (2p-l+6-a) 



e i+a+j,i+a ® e i+j,i- (22) 



The matrices and R~{z) are different; nonetheless they share many properties, 

including that they square to the identity and obey the unitarity condition. Moreover, 
we have the following proposition: 

Proposition 5.1. The descendant R~ (z) of Equation (22) satisfies the Yang-Baxter 
equation if and only if R + (z) given by Equation (21) does. 



" ' z _|_ w (2p-l+b-a) 



Proof. We use the identity 

c+f 

n 1 1 + zw^ p - i+b -^ 

p=c+l v 

This implies that 

m— 1 

= (-i)tI V 



(-1) 



j,j,a,ft=0 



6j 



n i 



P =i 



z _|_ .y/sp-i+fc-a) 



We now consider a basis transformation. Given any A G C, there exists a basis 
transformation which yields 

We choose A such that our operator becomes (after scaling) 



1 



m— 1 



R-{z) = - V 



(6+f)j 



n i 



P =i 



z _|_ w (2p-l+6-a) 



i,j,a,b=0 

We define the functions 

m— 1 a 

9(a,j)(z) = 

b=0 p=l 

which allow to write the operators 

m— 1 

R + {z)= ^ g(a,j)(z)e i+a+ j ti+a ® e i+ j ; 



, r -I- ■ ? ,,( 2 P~ 1 + b ~ a ) 

V^TT - 

11 \ 1 _|_ 2| y(2p-l+b-a) I ' 



i,j,a=0 



and 



m— 1 



i,j,a=0 
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Hence R~(z) differs from R + (z) only by a basis transformation and a permutation of 
the entries. We calculated previously that R + (z) satisfies the Yang-Baxter equation 
if and only if 



m—l 



[9(a,k-d)(x)g(k,a-b)(xy)g(b,c-k)(y) - g( c ,k-b)(x)g(k,c-d)(xy)g(d,a-k)(y)] = 0, 



k=0 



< a, b, c, d < m—l. Similarly we have that R (z) satisfies the Yang-Baxter equation 
if and only if 

m— 1 

X! [9(a+f,k-d){x)g(k+^.a-b)(xy)9(b+^,c-k)(y) - 9( c +!fM-b){x)g(k+!f.c-d){xy)g(d+f,a-k){y) = 0, 

fe=0 

< a, b, c, d < m — 1. As we can consider each of the indices of the functions modulo 
m we find that the two above constraints are equivalent, hence the result. 

□ 

Thus we have found another family of conjectured descendants. We explain the 
existence of R~ (z) by considering the values of functions associated with the irreps. To 
obtain R + (z) we simply set every function associated with irreps from the conjugacy 
class {e} equal to one, i.e. /( ,b) = 1 for b E Z. Conversely, R~(z) is obtained by 
setting every function associated with irreps from the conjugacy class {a m } equal to 
one, i.e. /(^^(l) = 1 for b e Z. We recall that e and a m are the central elements in 
D 2m when m is even. It is in part due to this enlarged centre that we obtain more 
general solutions in this case. 

It is also possible for us to ignore the unitary condition and only impose that R(z) 
is self-adjoint. This leads to the constraints 

f(o,2b)(z) = f(o,o)(z) = (f(o,o)(z))* and f(o, 2 b+i)(z) = f(o,i)(z) = (/(o,i) (*))*, 
for all b E Z and zel. Using these constraints we set 

f(o,o)(z) = 1 + f(z) and /( ,i)(*) = 1 - f(z), 
where f(z) is an arbitrary real function. This yields the operator 

R{z) = R + {z) + f{z)R-{z). (23) 

This is invertible provided 

/(*) ± ±1- 

The function f(z) is equivalent to a second parameter. That is, the operator 

R(z,fj) = R + (z)+ fiR-(z), (24) 
is invertible for /j, ^ ±1. Furthermore, we have the following: 
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Proposition 5.2. The operator R(z , ji) given by Equation (24) satisfies 

R 12 (x, X)R 23 (xy, n)Ri 2 {y, v) = R 23 (y, v)R l2 (xy, n)R 23 (x, A), (25) 

if and only if R(z) given by Equation (23) is a solution to the Yang-Baxter equation 
with no constraints on f(z). 

Proof. Consider a solution of the Yang-Baxter equation R(z) G End(V n <S>V n ) which 
contains an arbitrary function f(z). Furthermore suppose every entry of R{z) can be 
written as a polynomial in terms of z and f(z). Let 

Q = Ri 2 (x)R 13 (xy)R 23 (y) - R 23 (y)Ri 3 (xy)R 12 (x). 

We observe that every entry of Q is expressible as 

oo 

E ^ k (x,y)f(x)f(xy)f(y), 

i,j,k=0 

where h\j k (x,y) are polynomials in x and y and I indexes the entry of Q. As R(z) is 
a solution of the YBE, i.e. SI = 0, and f(z) is an arbitrary function, we deduce that. 

h ijk( x ,y) = o, V x >y e c /{°> s -t- x + y- 

Let y E C/{0, 1, oo}. There are at most we have four values of x for which h l ijk (x, y ) 
can be non-zero, but h\j k (x,yo) is continuous in x. This means that 

h[ jk (x,y o )=0 VxeC. 

By symmetry 

h\ jk (x, y) = 0, V(x, y) G (C x C)/({0, 1, oo} x {0, 1, oo}). 

Thus there are at most 9 points in which h\j k (x, y) can be non-zero; however h\j k (x, y) 
is continuous in x and y. Hence 

h ijk( x ,y) =°, Vx,yeC. 

Let R(z,fx) be the operator derived from R(z) in which we have replaced the free 
function f(z) with fj,. Every entry of R(z) must be expressible as a polynomial in 
terms of z and \i. If we let 

Q = R 12 {x, X)R 13 (xy, n)R 23 {y, v) - R 23 (y, v)Ri 3 (xy, n)R 12 (x, A), 

x, y, A, /i,i> e C, then we find every entry of f2 can be written 

oo 
i,j,k=0 
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As shown previously 

h ijk( x >y) = 0, Vi,|/GC. 
Thus f2 = and R(z, //) satisfies 

i?i 2 (x, X)R 13 (xy, n)R 23 (y, v) = R 23 (y, v)Ri3A x v) R M x , 

Conversely, suppose we have /i) which satisfies 

R 12 (x, X)R 13 (xy, fi)Rv3(y, v) = R 2 3(y, v)Ri3(xy, n)Rvi{x, A). 

If we consider R(z) = R(z, f(z)) where f(z) is an arbitrary function then R(z) must 
satisfy the YBE. To recover the result, we use the fact that R(z) = P R(z) and 
R(z,fj) = P R(z,fj). This proof works if the entries of R(z) can be written as a 
polynomial f(z) whose coefficients are rational functions of z. We are able to scale 
R(z) by the product of all the denominators of the coefficients of f(z), hence turning 
it into a polynomial. 

□ 

We are able to verify using Maple that R(z) given by Equation (23) satisfies the 
Yang-Baxter equation for even n up to 12, and hence that R(z, /i) satisfies Equation 
(25). We conjecture that this holds true for all even n, and have shown it to hold in 
the limit x = y = 0. 

6 Connection to the Fateev— Zamolodchikov model 

Closely associated with the Yang-Baxter equation is the star-triangle relation (STR), 
given by 

N-l 

W{x\a - d)W{xy\d, c)W{y\d - b) = W{x\b - c)W(xy\a - b)W{y\a - c), 

d=0 

for < a, b, c < N — 1. One well-known solution of the STR is the iV-state Fateev- 
Zamolodchikov model [23], which has weights 

3=1 ' 3=1 

for < I < N — 1, where A is a primitive 2iVth root of unity. These weights are 
extended by the relations 

W(z\l) = W(z\N + 1) and W(z\l) = W(z\N + I), 

while also satisfying 

W(z\l) = W(z\N -I) and W{z\l) = W{z\N - I), 
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for I e Z. These weights satisfy the STR and lead to the i?-matrix defined by 



TV 

R(x,y) = Rb a?a 2 {x,v)e b 1,01 eb 2) a 2 , 

01,02,61,62=1 



where 



and 

= ^(^rVi - fc 2 )W(^yr 1 |a 1 - ^(r^h - - 60- 

Through a private communication with V. Bazhanov and J. Perk [30] we learnt 
of a connection between the D(D n ) solution and the Fateev-Zamolodchikov model. 
Specifically in the private communication a limiting case of the 3-state Fateev- 
Zamolodchikov model was shown to reduce to the D{D^) solution. Using the ideas 
presented in [30] we are able to establish a connection between the TV-state Fateev- 
Zamolodchikov model and D(D n ) (or equivalently D(D 2n )) solution in the case where 
N = n and n is odd. 

To investigate which limit of the Fateev-Zamolodchikov model leads to the D(D n ) 
model, we determine when the Fateev-Zamolodchikov i?-matrix squares to the iden- 
tity. Using certain properties of the weights it is possible to show that 

R(x,y)R(x~ T ,y~ T ) cx I <g> 7, 

where 



x = I Xl ) and x T = [ X \ 
\x 2 J \ Xl 

Thus the inverse of the i?-matrix is known up to a scalar multiple. Furthermore we 
are able to show that if the i?-matrix squares to the identity then it is equivalent to 
one which satisfies the constraint 

x = x~ T and y = y~ T . 

This implies that the i?-matrix can be reduced to 

N 

y)= ^oit* ( x > y) e bi, ai ® e b2 ,a 2 , 

01,02,61,62=1 

where 

RaX&y) = W{xy- l \ ai - b 2 )W(x~ 1 y- 1 \a 1 - a 2 )W(x- 1 y\a 2 - &i)W(a;# 2 - 61). 
To obtain the difference property we set 

R(z) = lim R(x,y), (26) 

x,y— >oo 
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where z — ^. For odd n < 11 we are able to computationally verify that this R- 
matrix is equivalent (up to a basis transformation) to the D(D n ) i?-matrix while 
setting A = — w~ x . For odd n > 11 we can verify that -R(O) is indeed equivalent to 

We now briefly comment on the descendants obtained from D(D 2n ), n even. There 
are two distinct _R-matrices, R + (z) and R~(z), which both square to the identity. 
The multiplicities of the eigenvalues of R + {z) and R~{z) differ; furthermore neither 
eigenvalue spectrum matches that of R{z) as defined by Equation (26). Hence unlike 
the case when n is odd, the i?-matrix (23) is not equivalent up to basis transformation 
of a limit of the Fateev-Zamolodchikov i?-matrix. 



7 Summary 

In this paper we used the framework of descendants to construct i?-matrices from the 
Drinfeld doubles of dihedral groups. For 3 < n < 17 and w a primitive nth root of 
unity, 



n— i 



R(*)= E 

i,j,a=0 L b=0 p=l 

is a descendant of the six-vertex model 



n—1 a / 



z _|_ ^(Sp-l+b-a) 



+ zw( 2p ~ 1+b ~ a } 



( wz 1 — w l z 



r{z) = 



\ 



o 






z 1 — z w — w 1 
w — w' 1 z~ x — z 



&i+j,i+a ® &i+a+j,i 









wz 



1 - w~ x z j 



with corresponding L-operator 



n-l 



L ( z ) = E K^ 1 ' 2 + w ie 2,i) ® e M + z [ei.i ® + e 2 , 2 ® e i+M ]} . 



j=0 



We conjecture that this holds true for all n > 2. 

We also showed that when n is even we obtain i?-matrices with a second, non- 
spectral parameter. Specifically, let n > 2 be an even integer and w a primitive nth 
root of unity. Given 



n-l 



R + (z) = E 

i,j,a=0 



n— 1 



n 2^ ^ 1 1 I 1 + 2 W (2p-l+6-a) 

{,=0 p=l v 



and 



n-l 

E 

i,j,a=0 



n-l 



V ^ 11 \ 1 + zw (2p-l+b-a) 



b=0 



p=l 



&i+j,i+a ® &i+a+j,ii 
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then 

R(z,fi) = R + (z) + fiR-(z) 

satisfies Equation (25) for n < 12 and is conjectured to satisfy it for larger n. More- 
over, R(z,fx) and L(z) together satisfy 

L 1 2(x)L 13 (y)R 2 3(x~ 1 y,n) = R 23 (x~ 1 y, fx) L 13 (y) L 12 (x) . 
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